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Abstract 

We present a general discussion of radiative four-meson processes to 0(p 4 ) in chiral pertur- 
bation theory. We propose a definition of "generalized bremsstrahlung" that takes full advan- 
tage of experimental information on the corresponding non-radiative process. We also derive 
general formulae for one-loop amplitudes which can be applied, for instance, to rj — > 3tij, 
mi — > TTTT^f and K — > 37r7. 
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1. Chiral perturbation theory (CHPT) |T|, |2j, ||| incorporates electromagnetic gauge invari- 
ance. To lowest order in the derivative expansion, 0(p 2 ) in the meson sector, amplitudes 
for radiative transitions are completely determined by the corresponding non-radiative am- 
plitudes. Direct emission, carrying genuinely new information, appears only at 0(p 4 ). Nev- 
ertheless, even in higher orders of the chiral expansion part of the amplitude is related to 
the respective non-radiative process. In order to isolate the direct emission amplitude, an 
operational definition of "generalized bremsstrahlung" is needed. 

In this letter, we investigate in a general manner radiative transitions involving four 
pseudoscalar mesons and one real photon. Possible applications to be discussed elsewhere 
include rj — > 3ir'~f and tttt — > 7T7T7 in the strong sector and the nonleptonic weak decays 
K — > 3ttj. Our purpose is twofold: 

i. We extend Low's theorem Q by terms of O(k) {k is the photon momentum) to define 
generalized bremsstrahlung. This part will include in particular all local terms of 0(p 4 ) 
that contribute also to the non-radiative four-meson transition. 

ii. We give a compact expression for the loop amplitude of a general four-meson process 
with a real photon. The resulting formula is immediately applicable to both strong and 
nonleptonic weak processes. We also consider the limiting case of a radiative three- 
meson amplitude to recover known results for K — > 27r7 decays || §, 0, ||. 

2. The amplitude for a four-meson transition with a single photon can be decomposed into 
an electric and a magnetic part: 

A{ip a ip b ip c ip dl ) = ee^fc)^ + e^ pa M vpa ) (1) 

with 

fc"£„ = 0, e^W" = . 

The magnetic amplitude M vpc can only occur in nonleptonic weak processes (for a review, see 
Ref. 0]). It appears first at 0(p 4 ) as a tree-level contribution and will not concern us further. 
Here, we are only interested in the electric amplitude that is in particular sensitive to 
bremsstrahlung. The kinematics of the process is specified by five scalar variables which we 
choose as 

s=(pi+P2) 2 , v = p 3 (p 1 -p 2 ), U = k-pi (z = l,...,4) (2) 

with 

4 

Y /Pi + k = 0, t x + t 2 + t 3 + u = o . 

i=l 

Any three of the U together with s and v form a set of independent variables. 

The non-radiative transition is characterized by the two Dalitz variables s and v. De- 
noting the non-radiative amplitude by A(s, u), Low's theorem M amounts to the following 
expansion in the photon momentum k: 

= A(s,v)W 

n dA(s,v) KU dA(s,v),. u . UN 
+ 2 A i2 + - (A13 - A23) 

+ 0(k) S (3) 
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with (the meson charges in units of e are denoted qA 



4 /J, 
i=l 

M', = A£ = (qitj - qjU) B% 

D?. = — £) - = £L - £l . (4) 

The explicit terms in (^) are often called "internal bremsstrahlung" . It is straightforward to 
show that there are no terms of 0(k) at lowest order in the chiral expansion. Thus, for ra- 
diative four-meson processes the leading chiral amplitude of 0{p 2 ) is completely determined 
by the non-radiative amplitude A(s, v) as expressed by Eq. (§). 



3. At 0(p 4 ), there are as usual both one-loop and tree-level contributions with a single 
vertex from the strong Lagrangian £4 0, |3|] or the nonleptonic weak Lagrangian C^ s=1 
[ [10| , [TTf| . Let us first consider the tree-level amplitude. The different terms in either £4 or 
Cf s=1 that can contribute to the processes under consideration can be grouped in four 
classes: 



A. Terms of 0(mV) without derivatives: fully covered by internal bremsstrahlung (|3|). 

B. Terms of 0(m q ) with two (covariant) derivatives: again included in (Q). 

C. Four-derivative terms: in general not fully covered by 

D. Terms with two derivatives and one field strength tensor containing the electromagnetic 
field: contribute only to the radiative transition and thus are never included in 



Obviously, groups A,B correspond to internal bremsstrahlung while the contributions 
of type D belong to direct emission. Class C falls in between if we adopt Eq. (||) as the 
definition of bremsstrahlung. On the other hand, it would have both conceptual and practical 
advantages to include all terms under the heading "bremsstrahlung" that contribute to both 
radiative and non-radiative transitions. One practical advantage arises for K — > 37r7 decays 
where the low-energy constants of the four-derivative terms are only partly known Jl2|, [TT| . If 
those terms could be included in bremsstrahlung, we may use experimental data for K — > 3n 
decays directly without having to worry about the values of the aforementioned coupling 
constants [13| . 

In order to incorporate class C in what we shall call generalized bremsstrahlung, we must 
add explicit terms of 0(k) to Low's formula ([3D. The clue for the solution is the observation 
that both £4 and C^ s=1 give rise to at most three independent four-derivative couplings at 
the mesonic level: 



D fl (p a D^(p b D u if c D u (p d , D fl (p a D' M !f c D u ip b D u (f d , D fl ip a DL l ip d D u ip b D u <p c 

DpWa = (<9 M + iqaeA^Lpa . (5) 
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At the same time, we have three independent second derivatives of the non-radiative am- 
plitude A(s, v). Therefore, the following extension of (|3]) solves the problem (GB stands for 
"generalized bremsstrahlung" ) : 

E» = E» B + 0(k) 
E» B = A( S ,^ + 2^^A^ + ^&^(A-3-A- 3 ) 

+ 2 ^^l {tl + h)K + l^j^iih - t 2 )(A? 3 - A£ 3 ) - htJT] 

+ 2 ^r foK- tl A% 3 ]. (6) 



It is important to realize that Eq B in (^) does not contain all terms of at most 0(k). 
In fact, it is impossible in general to relate all terms of O(k) to derivatives of A(s, v). On 
the other hand, the definition of generalized bremsstrahlung in (J6|) guarantees that all local 
(counter)terms that contribute to both radiative and non-radiative processes (classes A,B,C) 
are included in Eq B . The difference E^ — Eq B is at least O(k) and will be referred to as the 
direct emission amplitude (of the electric type). 



4. We now turn to the loop amplitude. Most of the renormalization procedure can trivially 
be carried over from the non-radiative to the radiative amplitude because all diagrams 
(tadpoles) relevant for mass, charge and wave function renormalization contribute only to 
internal bremsstrahlung. Thus, this part is completely taken care of by the non-radiative 
amplitude. Moreover, for a real photon there are no diagrams of the form factor type where 
the photon emerges from a mesonic bubble. 

Restricting attention to transitions of at most first order in the Fermi coupling constant, 
the only non-trivial diagram is of the type shown in Fig. 1 where the photon can hook on to 
any charged meson line and to any vertex with at least two charged fields. The two vertices 
are either both from the lowest-order chiral Lagrangian C 2 0, |3j (strong transition) or one 
from £2 and one from C^ s=1 ( c f-> e -§-> R- e f- 10) f° r a nonleptonic weak transition. Despite 
the comparative simplicity, diagrams of the type displayed in Fig. 1 with a photon in all 
possible places generate a considerable number of terms due to the derivative structure of 
vertices. Moreover, there are usually several permutations (1234) — > (abed) that have to be 
added for a given process. 

We have therefore found it useful, both for our own work and for possible future applica- 
tions, to present the loop amplitude in a compact form suitable for any strong or nonleptonic 
weak transition. For this purpose, we first calculate the loop contribution to the non-radiative 
amplitude A(s, v). 

We characterize the vertices Vi, V 2 in momentum space by real constants a^, bf. 

Vi = a + a lPa - Pb + a 2 p a -x + a 3 (x 2 - Ml) + a 4 (y 2 - M 2 y ) + a 5 (p 2 a - M 2 ) + a 6 (p 2 b - M 2 ) 
V 2 = bo + b^-p. + hp.-x + b^-MD + b^-M^ + hipl-MD + bM-M 2 ) . 

(7) 
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Figure 1: One-loop diagram for the four-meson transition. For the radiative amplitude, the 
photon must be appended to every charged meson line and to every vertex with at least two 
charged fields. The vertices Vi, V 2 are defined in Eq. ([?]). 



In calculating the loop amplitude for the non-radiative process, we do not associate the vari- 
ous scalar products with the Dalitz variables s, v that will depend on the specific assignment 
(abed) — f (1234). Instead, we write the loop amplitude as a function of the four-momentum 

P = Pc + Pd • (8) 

The non-photonic loop amplitude of Fig. 1 can be represented in the following form that 
will turn out to be useful (all external lines are on-shell): 

F(P) = A(M X ) [a^pa-p,, + a 4 b lPc -p d + a A b 4 (P 2 + M 2 - M 2 ) + a 6 4 + a 4 6 ] 
+ A(M y )[a 1 b 3 p a -p b + a 2 b 3 p a -P + a 3 b x p c -p d + a 3 b 2 p c -P 

+a 3 b 3 (P 2 ~ M 2 + M 2 ) + a b 3 + a 3 b ] 
+ B(P 2 , M x , M y )[a b + a bip c -p d + aib p a -p b + a x bxp a -PbPc-Pd\ 
+ B 1 (P 2 ,M x ,M y ) [a b 2 p c ■ P + a 2 b p a • P + a x b 2 p a ■ p b p c ■ P + a 2 b x p c ■ p d p a ■ P] 
+ a 2 b 2 [p a -p c B 20 (P 2 ,M x ,M y )+p a -Pp c -PB 22 (P 2 ,M x ,M y )] . (9) 

The various functions in @ are as defined conventionally (in d dimensions): 

1 r d d x 

A(M) -- 



i J (2n) d x 2 — M 2 

(B,B x P^ gilv B 2 , + P,P v B 22 ) = 7 / (t^ ( X 2 - Ml)[{x -P) 2 - M 2 ] ' (10) 
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Several comments are in order at this point. 



i. In order to limit the number of terms generated, it is preferable to express the functions 
Bi, B 20 and B 22 via the usual recursion relations in terms of A, B only at the very end. 

ii. The analytically non-trivial part of (0), involving the various B functions, contains 
only the on-shell couplings a ,a>i,a 2 ,b ,b 1 ,b 2 . The off-shell couplings a 3 , 04,63,64 ap- 
pear only together with the divergent constants A(M). Since these terms are polyno- 
mials in the momenta of at most degree two, they will enter in the radiative amplitude 
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only through internal bremsstrahlung. All the divergences in will be absorbed by 
counterterms belonging to classes A,B,C of the previous classification. As emphasized 
before, the generalized bremsstrahlung contains all these divergences plus the cor- 
responding counterterms. 

iii. For a reason soon to become evident, we have chosen to express F(P) in terms of the 
scalar products 

Pa-Pb, Pc-Pd, P 2 , Pa-P, Pc-P, Pa-Pc ■ (H) 

In other words, we have not used kinematical relations to write F(P) in terms of only 
two independent scalar variables [like s, v defined in Eq. (§)]. 



5. We now have all the ingredients for calculating the radiative loop amplitude E^ cor- 
responding to the diagrams of Fig. 1 with a photon in all possible places. For the general 
vertices Vi,V 2 given in (]?]), this loop amplitude contains several hundred terms even be- 
fore reducing the various B functions via recursion relations. A compact representation will 
therefore be of great use for avoiding tedious repetitions of the same procedure. 

We find it useful to decompose the radiative loop amplitude into two parts: 



-^loop 



(12) 



The more tedious part of the calculation is contained in the amplitude that can be 
expressed through various derivatives of the non-radiative loop amplitude F in (§) with 
respect to the scalar products (|TT|). I n some of the following terms, the momentum P has to 
be replaced by P + k, leaving all scalar products unchanged that do not contain P explicitly: 



F{P)?y - 
dF 



F(P + k) 



+ 



d{ Pa -P) 



k-P 



(P)Kp + 



dF 



d(p a -Pb) 



dF 



d(p c -Pd) 



(P + k)^ cd + 



(P)K b 

dF 
d(p c -P) 



(P + fc)A£ 



cP 



q a t c 



dF 



d(p a -Pc 



-{P)-q c t a 



dF 



d(p a -Pc 



■{P + k) 



ac 



- i^sic + qd)t a t c 



d 2 F 



d(p a -P)d( Pc -P) 



d 2 F 



aP 



d(p a -P)d( Pc -P) 



(13) 



We have used the definitions @. When P appears as an index (e.g., in A^ P or D^p), the 
corresponding momentum and charge in (H) are P and q c + Qd, respectively. For better 
understanding of the notation in (13), we give two explicit examples: 

dF 



d(p c -P) 
d 2 F 



d( Pa -P)d( Pc -P) 



(P + k) = a 3 b 2 A(M y ) + (a b 2 + a 1 b 2 p a -p b )B 1 ((P + ky,M X} M y 
+ a 2 b 2Pa -(P + k)B 22 ((P + k) 2 , Mb, M y ) 
(P) = a 2 b 2 B 22 (P 2 ,M x ,M y ) . 



(14) 
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The second part if M of the loop amplitude (|12[ ) cannot be expressed in terms of F or 
derivatives thereof. Since the dominant contributions to Ey oop are usually due to pion loops 
(if they contribute at all), we give the explicit expression for only for the case of equal 
loop masses (M x = M y =: M). In this special case, if M takes on the following compact form: 

H» = a 2 (t b p^-t a p^{(q x -q^(2b + 2b lPc -p d + b^ c -P)C w (P^-k-P)^ 

+ h(q x + qy^pc-P^P 2 , -k-P) + 2t c Cl 2 {P\ -k-P) - p c -PC^ (P 2 , -k-P)]} 
+ &a(*dP£ - tcP%){(q x - q y )[2a + 2a lPa -p b + a 2 (p a -P + t a )]C^{{P + k) 2 , k-P) 
+ a 2 (q x + q y ) [-2{p a -P + t a )C7i((P + k) 2 , k-P)- 2t a Gf 2 ((P + k) 2 , k-P) 
- (Pa-P + t a )C^({P + k) 2 , k-P)]} . (15) 

The functions Cy are defined as 

Q,(u,v) = C - M - C ^ 0) (16) 

v 

in terms of the three-propagator one-loop functions CV,(p 2 , k-p) for k 2 = 0: 

IP /™fd ,-jsi f />-» />■ rr* r f' r T' ~\. 

77 (2vr) d (x 2 - M 2 )[(x + p) 2 - M 2 ][{x + k) 2 - M 2 ] = 
= {C 20 (p 2 , k-p)g llu + . . . , C 3 i(p 2 , k-p){p p g vp + p v g w + p p g pv ) 

+ C 32 (p 2 , k-p){k p g vp + k u g w + k p g pu ) + ...}. (17) 

As in the case of F(P) in @, it is advisable not to use the standard recursion relations for 
the functions C 2 o, C31 and C32 in ([T5) until the actual numerical analysis. At the expense 



of introducing the functions B, B 20 , B 22 defined in (fLOf), one may express the functions 
Cij((P + k) 2 , k-P) in terms of the Cij(P 2 , —k-P) or vice versa. 

The following comments (valid also in the case of different loop masses) explain the 
motivation for splitting the loop amplitude Ej 1 in flT2| ) into two parts. 

i. The amplitudes in (|T3| ) and if M in (plj) are separately gauge invariant. 

ii. The amplitude is finite and at least of 0(k) as is evident from Eqs. fll5|) , (pT6|) and 
0T7|). Moreover, it only contains the on-shell couplings ao, a%, a 2 , bo, b\, b 2 defined in (|7|) 
and the charges q x , q y of the particles in the loop. Of course, we have 

Qx + Q y = -q a ~Qb = Qc + Qd ■ (18) 



iii. The amplitude contains the generalized bremsstrahlung (^j) for the non-radiative 
loop amplitude @. If we denote by -Eq B (1oop) the result obtained by inserting for 
A(s, v) the on-shell loop amplitude (|9|) in Eq. (||), then the difference 

A" = G» - ^ B (loo P ) (19) 

is at least of O(k). Moreover, by construction of Eq B the divergences in A M are renor- 
malized by counterterms of class D only, i.e. by counterterms with an explicit field 
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strength tensor. In the strong sector, the relevant couplings of 0(p 4 ) are Lg for chiral 
SU(3) [£ e for SU(2)] and N u , Ni 5 , iV 16 , N 17 for the octet part of the nonleptonic weak 
Lagrangian [|TTJ . Finally, it can be shown that if CI2&2 = then A^ is finite and at least 
of 0(k 2 ) for s = £ 4 =1 Mf/3, v = and arbitrary U. 

iv. The apparent asymmetry of G M and under interchanges a <-> b or c <-> d is due 
to the asymmetric definition of vertices (|7]). For the same reason, C 1 and if M are in 
general not invariant under interchanges of the loop particles x and y. 

For the realistic case with experimental information on the non-radiative amplitude 
A(s, u), the complete electric amplitude to 0(p 4 ) accuracy can be written as 

E» = E£ B (exp) + < unter + ( A M + H») , (20) 

ioops 

where several loop diagrams may have to be added for a given physical transition. Only coun- 
terterms with an explicit field strength tensor must be included in £^unter • Consequently, only 
the renormalized coupling constants £g(^g) and/or iV[ 4 , . . . , iV[ 7 appear in (p0|). Of course, 
the amplitude is finite and scale independent by construction. All other counterterms of 
0(p A ) are hidden in i?Q B (exp). An alternative approach is to use instead of the experimental 
amplitude i?Q B (exp) in (ED|) the theoretical prediction Eq-q (theory) in terms of the amplitude 
A(s, v) calculated in CHPT to 0(p 4 ) accuracy. Both approaches are equivalent to 0{p A ). The 
difference between them gives an indication of the size of effects of 0(p 6 ) and higher. 

6. As a final application, we consider the limit q a — > 0, p a — ► to connect with known 
results for K — ► 27r7 decays. In this case, the vertex V\ is necessarily of the weak nonleptonic 
type and the coupling constants 01,02,05 disappear. A straightforward calculation shows 
that G* 1 becomes in this limit 

G" = F(P + k)?T (21) 

where F(P + k) = F(—pb) is now the on-shell loop amplitude for the decay b — > c + d. 
Eq. ([H]) is nothing but the familiar bremsstrahlung amplitude for a radiative three-meson 
transition and as a consequence 

A^ = . (22) 
Likewise, the amplitude in ( |T5D reduces in the three-meson limit to the single term 

H» = 2a b 2 (q x - q y )C 20 (Ml -k- Pb ){t dP ^ c - t c f d ). (23) 

Thus, in this limit, the loop contribution to the direct emission is finite and proportional to 
a , the on-shell tree-level amplitude for the nonleptonic weak transition b —>■ c + d. 

7. To summarize, we have presented a general discussion of radiative four-meson processes 
to 0(p 4 ) in CHPT. We have proposed a definition of generalized bremsstrahlung in Eq. 
that has the advantage of including all counterterms of 0(p 4 ) that contribute to both radia- 
tive and non-radiative amplitudes. For general vertices of 0(p 2 ) that encompass all strong 
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and nonleptonic weak transitions of interest, we have calculated the non-trivial loop ampli- 
tudes in terms of two gauge invariant parts. The amplitude G M given in (|T3|) is expressed in 
terms of the non-radiative loop amplitude F in @. The divergences are all contained in G 11 . 
The remainder if M given in Eq. fll5|) is finite and at least of O(k). 

Applications of these general results to K — > 37r7 decays and other radiative four-meson 
processes will be presented elsewhere [ |I~3|j . 
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